Resonant plasmonic effects in periodic graphene antidot arrays 
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We show that a graphene sheet perforated with micro- or nano-size antidots have prominent 
absorption resonances in the microwave and terahertz regions. These resonances correspond to 
surface plasmons of a continuous sheet "perturbed" by a lattice. They are excited in different 
diffraction orders, in contrast to cavity surface plasmon modes existing in disconnected graphene 
structures. The resonant absorption by the antidot array can essentially exceed the absorption by a 
continuous graphene sheet, even for high antidot diameter-to-period aspect ratios. Surface plasmon- 
enhanced absorption and suppressed transmission is more efficient for higher relaxation times of the 
charge carriers. 

PACS numbers: 42.25.Bs, 41.20.Jb, 42.79.Ag, 78.66. Bz 



Graphene now is widely known not only for its 
amazing transport properties^ but also for a signif- 
icant potential in photonics P Most applications re- 
quire that the electromagnetic (EM) field is first ab- 
sorbed by graphene. In a large part of the EM spec- 
trum the two-dimensional conductance of one mono- 
layer a is small, so that a = 2na /c <C 1. Then the 
absorption coefficient A for a free-standing graphene 
can be approximately written at normal incidence as 
A ~ 2Re(a). In the optical regime graphene behaves 
as an absorbing dielectric (its effective refractive index 
is a purely imaginary number): a opt = 7rao/2, with 
ao = 1/137 being the fine-structure constant. This 
leads to the classical result A opt = nao ~ 2.3%. In 
both terahertz (THz) and microwave frequency ranges 
the situation is rather different. The conductivity 
is contributed predominantly by the intra-band elec- 
tronic transitions, so that graphene behaves like a 
two-dimensional metal. For low temperatures its con- 
ductance can be approximated by a Drude-Lorentz 
modeP^axHz = 2iao\/J>\/h(uj + ir -1 ) depending upon 
the chemical potential /i, frequency uj and relaxation 
time r. In this case for ujt <C 1 the absorption coeffi- 
cient reads Athz — 4ao|/i|/(^™ 2 ) and can become of 
order of tenths of a percent. Obviously, it is challeng- 
ing to use graphene as an efficient absorber. 

Fortunately, one of the promising points of graphene 
is that it supports surface waves [graphene surface 
plasmons (GSPs)].™£J GSPs allow for a strong sub- 
wavelength confinement of the EM fields and thus 
for a much higher absorption. Previously it has 
been demonstrated that periodic arrays of graphene 
ribbons)- 12 * 14 * discs^ and rings^can resonantly absorb 
THz radiation due to excitation of GSP cavity modes. 
However, this kind of structures is not "electrically- 
continuous" preventing therefore a direct current pass- 
ing through the samples. An alternative solution can 
be a graphene periodic antidot array (GPAA)P^31 
which presents both continuity (to provide the electric 



transport) and periodicity (to transfer the momentum 
to GSPs from the incident wave). The antidots in a 
graphene sheet can be considered as regions of zero 
conductivity due to either (i) actual holes in graphene, 
(ii) the influence of an external voltage or (iii) appro- 
priate chemical doping. 

GPAAs could open interesting possibilities for elec- 
tromagnetic control via GSP-assisted effects. This can 
be an important point for designing photo-electric de- 
vices, e.g. tunable absorbers that can act at the same 
time as conducting electrodes. 




FIG. 1: (Color online) The schematic of the studied sys- 
tem: a plane electromagnetic wave impinging onto a peri- 
odic antidot array in a graphene monolayer placed onto a 
dielectric substrate. 

In this Letter we report on GSP-enhanced ab- 
sorption and suppressed transmission in continuous 
GPAAs with circular antidots. We characterize reso- 
nances, presenting various dependencies of absorption 
A, transmission T, and reflection R coefficients upon 
several parameters: period L, antidot diameter d and 
relaxation time r. We consider a square GPAA placed 
on the substrate with dielectric permittivity e. The 
system is illuminated by a normal-incident monochro- 
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matic plane wave with wavevector k from the vacuum 
half-space. The electric field of the incident wave is 
parallel to one of the Bragg's vectors of the lattice. The 
scheme of the structure is shown in Fig. [I] Graphene 
is modeled by a conductivity cr, computed within the 
random-phase approximation,^"^ which is valid pro- 
vided that the mean free path is shorter than the su- 
perlattice scale. In this work we assume room temper- 
ature, T = 300 K. 
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FIG. 2: (Color online) The spectra of absorption, reflection 
and transmission coefficients for square GPAAs. (a) A, R, 
T for the suspended GPAA and GPAA on the substrate 
with e — 3 and r = 40 ps. Rf and Tp corresponding to the 
continuous graphene sheet (Fresnel coefficients) are shown 
by the discontinuous curves, (b) Absorption spectra for 
GPAAs at different r for s = 1. The spectra for continuous 
graphene monolayers are shown by the dashed lines. The 
insets to (b) show the colorplots for the spatial distribution 
of the electric field: modulus (below) and the real part of 
the projection on the direction of the incident wave electric 
field shown by an arrow (above). In all panels \i — 0.2eV, 
L — 2d — 5/xm. 

An example of the wavelength spectra for A, T and 
R is shown in Fig. [2] (a) corresponding to L = 2d = 
5 jam. The simulations have been conducted by using 
the finite elements method (FEM)P^We choose this 
value of L for demonstration purposes, in order to have 
prominent resonances in the window 20 — 100/im; the 
dependence with L will be addressed later on. Both 



the case of free-standing GPAA and GPAA on the sub- 
strate are shown. While absorption and reflection co- 
efficients have a set of resonance maxima, the trans- 
mission spectra possess pronounced minima. The po- 
sitions of the resonances correspond approximately to 
the condition of GSP excitation in different diffraction 
orders, i.e. G\Jn 2 + m 2 = fc p , where G = 2tt/L is the 
shortest Bragg's vector, and k p is the GSP wavevec- 
tor. For convenience, these resonances will be further 
labeled as [n, m]. In suspended GPAAs with suffi- 
ciently small antidots a good approximation for k p can 
be k p ~ k\[\ — a -2 , with a = a • [1 — f (f ) 2 ] being 
the average value of the normalized graphene conduc- 
tance in PGAA. The characterization of the resonance 
positions for large antidots will be presented below. 
For graphene in a nonsymmetric dielectric environ- 
ment the dispersion relation reads — + — + 2a = 0, 

1 q z i q Z 2 

where ^1,2 = ^1,2 — k 2 /k 2 and ei^ are the dielec- 
tric constants of the surrounding media. 9 This equa- 
tion cannot be resolved in algebraic functions and 
therefore is treated numerically. The positions of the 
three highest resonance wavelengths estimated from 
the above "undressed" dispersion law are shown by 
the discontinuous vertical lines in Fig. [2] (a). The min- 
imum/maximum corresponding to the longest wave- 
length resonance [1, 0] is redshifted with respect to the 
unperturbed resonance wavelength, while the next two 
(lower-wavelength) ones, i.e. [1,1] and [2,0] are blue- 
shifted. In a symmetric configuration the maximal lim- 
iting value for A is 1/2, but it can be increase d, for 
instance, placing a metal layer below the GPAAfHUU 

In order to better visualize the resonance picture, 
the inset to Fig. [2] (b) shows the spatial field distri- 
bution for both electric field modulus and the real 
part of the electric field component along the polar- 
ization of the incident wave. The colorplots are taken 
at the wavelengths corresponding to [1, 0] and [1, 1] res- 
onances at A = 53.26 fim and A = 42.06 /im, respec- 
tively. Let us point out that the plasmonic wavelengths 
X p corresponding to a continuous monolayer for these 
values of A are X p = 6.61 fim and X p = 4.12 /im. 

Notice that GSP resonances in PGAA are reminis- 
cent to short-range surface plasmon resonances studied 
in hole arrays in ultra-thin metallic films ™ However, 
in case of metallic structures, absorption and reflection 
peaks (transmission dips) appear in the optical region. 

Owing to much higher values of rj 27 * 28 l suspended 
samples can present a special interest. Nowadays, plac- 
ing graphene on a substrate severely reduces the mo- 
bility. For this reason, in this Letter we mainly con- 
centrate on the free-standing GPAA. However, for aca- 
demic purposes, in order to see other effects of the sub- 
strate not related to the change of mobility, in Fig. [2] 
(a) we compare the spectra for symmetric and nonsym- 
metric configurations at the same value of r. From this 
comparison it follows that the spectra are similar, and 
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the only significant difference between them consists 
in a displacement of the resonance for GPAA to longer 
wavelength compared to the free-standing case. This is 
related to the higher values of the GSP wavevector for 
a nonsymmetric dielectric environment. Let us turn to 
the behavior of the absorption coefficient with respect 
to the change of r. 
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FIG. 3: (Color online) (c) Absorption spectra for GPAAs 
at different antidot sizes d. The main figure is for L = 5 /mi, 
while the inset is for L = 1/xm. The curve corresponding 
to d/L = renders A for the continuous monolayer, e = 3, 
r = 40 ps, ii = 0.2eV. 

In Fig. [2] (b) the absorption spectra for different re- 
laxation rates r are compared. It is seen that higher 
r provide higher visibility of the resonances, especially 
for lower wavelengths. An important message here is 
that the GSP-induced absorption in GPAAs can not 
only be higher than the absorption corresponding to 
lower r, but can also largely exceed the absorption in 
the continuous graphene monolayer. Since resonances 
for different r show the same behavior in function of 
geometric parameters of GPAAs, we will further con- 
centrate on the case of the maximal predicted value for 
r in continuous films 

The resonances are very sensitive to the geometrical 
aspect ratio d/L. In Fig.plthe evolution of the longest- 
wavelength resonance [1,0] with respect to the antidot 
diameter change is illustrated for L = 5 jam. When 
the antidot size increases, the resonance suffers a red- 
shift and the intensity of the absorption peak increases. 
However, in the interval of d/L = (0.3,0.8) the reso- 
nance position is quite stable. Further increasing d/L 
up to the values d/L ^ 1 leads to the fast displacement 
of the peak. A point to note is that even for large d/L, 
for which the graphene area is smaller than the one 
of the antidot one, the maximal value of the resonance 
peak remains hight. The same behavior takes place for 
other values of the period. As an example, the case of 
L = 1 /am is shown in the inset to Fig. [3j Notice that 
in contrast to the arrays of graphene ribbons or discs, 
the resonance position for GPAA at a fixed value of d 



is strongly dependent upon the period due to Bragg's 
origin of the GSP resonance. 
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FIG. 4: (Color online) Characterization of the resonances 
in GPAAs with respect to the aspect ratio d/L. (a) Maxi- 
mal value of absorption A max as a function of d/L at a fixed 
L for the resonances [1,0] and [1,1]. (b) The positions of 
the maxima for [1,0] and [1, 1] resonances at different d/L 
for a fixed L. (c) The positions of the maxima for [1,0] 
and [1, 1] resonances at different d/L for a fixed d. In (b), 
(c) continuous lines represent the estimated position of the 
resonance, while the symbols correspond to the full calcu- 
lations. In all panels e = 1, /a = 0.2 eV and r = 40 ps. 

In order to gain a further insight into the proper- 
ties of GSP resonances in GPAAs, let us analyze in 
more detail the dependencies of the absorption max- 
ima A max and resonance wavelength A max upon d/L. 
We consider GPAAs with the L values down to 1/am 
and values of d down to lOOnm. In Fig. [4] (a), (b) A max 
and A max both for [1, 0] and [1, 1] resonances are shown 
as a function of d, at several fixed periods. For the 
[1,0] resonance, A max quickly reaches high values, af- 
ter which it saturates, remaining practically constant 
for all periods. In contrast, A max for the [1,1] reso- 
nance has an optimum close to d/L = 0.7, with this 
value being independent upon L. 
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The dependency of A max upon d and L changes 
in different regions of d/L. When d/L < 0.1, the 
dispersion relation of the GP and therefore A max , is 
governed by the average value of a, see the coinci- 
dence of the points and the curves in the region of 
small d/L in Fig. [4] (b,c). However, for larger d/L 
there is a modification of the GSP dispersion branches. 
The degenerated resonance diffraction orders become 
strongly coupled (directly and also due to multiple 
scattering processes through both evanescent and ra- 
diative diffraction orders), leading to an increase of 
the bandgap width and a shift of the bandgap center. 
In other words, the dependence of A max deviates from 
the estimation based on the average conductance, see 
Fig.0(b,c). 

To summarize, in this Letter we have shown that 
graphene periodic antidot arrays (GPAAs) can pro- 
vide a strong electromagnetic response in both mi- 
crowave and THz regions, where GSPs are excited. 
The absorption peaks and transmission dips in the an- 
tidot structure are more pronounced for less absorbing 
graphene samples, characterized by higher relaxation 



rates of charge carriers. The system discussed here 
can be realized in CVD graphene, where large sam- 
ples can be routinely fabricated. 29 The resonances in 
antidot arrays are due to excitation of dressed GSPs 
under the Bragg's condition. This makes the difference 
with electrically-isolated graphene periodic structures 
(stripes, discs, etc.), where the resonances are due to 
cavity GSP-modes. The continuous graphene struc- 
tures possess thus a supplementary tunability parame- 
ter L. Additionally, GPAAs can have applications for 
devices, in which the same graphene layer can be both 
an electric conductor and radiation absorber. 
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